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ABSTRACT 

We investigate the nature of 4 young and low-mass open clusters (OCs) located in 
the 2 nd and 3 rd quadrants with near-IR 2MASS photometry (errors < 0.1 mag). After 
field decontamination, the colour-magnitude diagrams (CMDs) display similar mor- 
phologies: a poorly-populated main sequence (MS) and a dominant fraction of pre- 
MS (PMS) stars somewhat affected by differential reddening. Pismis5, vdB80 and 
BDSB96 have MS ages within 5 ± 4Myr, while the MS of NGC 1931 is 10 ± 3Myr 
old. However, non-instantaneous star formation is implied by the wider (~ 20 Myr) 
PMS age spread. The cluster masses derived from MS+PMS stars are low, within 
~ 60 — 180 Mq, with mass functions (MFs) significantly flatter than Salpeter's ini- 
tial mass function (IMF). Distances from the Sun are within 1.0 — 2.4 kpc, and the 
visual absorptions are in the range Ay = 1.0 — 2.0. From the stellar radial density 
profiles (RDPs), we find that they are small (R c < 0.48 pc, -Rrdp ^ 5.8 pc), especially 
Pismis5 with R c xt 0.2 pc and Rrdp ~ 1.8 pc. Except for the irregular and cuspy 
inner regions of NGC 1931 and Pismis5, the stellar RDPs follow a King-like profile. 
At ~ 10 Myr, central cusps - which in old clusters appear to be related to advanced 
dynamical evolution - are probably associated with a star- formation and/or molecular 
cloud fragmentation effect. Despite the flat MFs, vdB 80 and BDSB 96 appear to be 
typical young, low-mass OCs. NGC 1931 and especially Pismis5, with irregular RDPs, 
low cluster mass and flat MFs, do not appear to be in dynamical equilibrium. Both 
may be evolving into OB associations and/or doomed to dissolution in a few 10 7 yr. 

Key words: ( Galaxy:) open clusters and associations: general; ( Galaxy:) open clus- 
ters and associations: individual: Pismis5, vdB80, NGC 1931 and BDSB 96. 



1 INTRODUCTION 

The first few ten Myrs represent the most critical phase in 
the l ife of a star cluster , to the point that only about 5% 
(e.g. lLada fc Lada| [2003) of the embedded clusters are able 
to dynamically evolve into bound open clusters (OCs). The 
rapid gas removal by supernovae and massive star winds 
associated with this period can produce important changes 
on the primordial gravitational potential. Obviously, this ef- 
fect depends essentially on the star-formation efficiency, the 
actual mass of primordial gas converted into stars and the 
mass of the more massive stars. 

Because of the rather rapidly-reduced potential, a sig- 
nificant fraction of the stars - the low mass ones in particular 
- moving faster than the scaled-down escape velocity may be 
driven into the field. Over a time-scale of 10 — 40 Myr, this 
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process in t urn, can dissolve most of the very young star 
clusters fe.g. iGoodwin fc~B astian 2006). 

Observationally, low-mass star clusters younger 
than about 10 Myr present an under-populated, de- 
veloping main sequence (MS) and a more conspic- 
uous population of pre-M S (PMS) stars. Typical 
exam ples are NG C 6611 (Bonatto. Santos Jr. fc Bical 
l2006h . NGC 475 5 (iBonatto et all l2006bl h" NGC 2239 



and NGC 2244 jBonatto fc Bicd l2009bh and Bochuml 
l|Bica. Bonatto fc Dutrall2008l '). On a large scale, the impor- 
tant changes to the potential that affect the early cluster 
spatial structure should also be reflected on the stellar 
radial density profile (RDP). Bochuml and NGC 2244, for 
instance, can be taken as representative of this scenario, in 
which an irregular RDP cannot be represented by a cluster- 
like (i.e. an approximately isothermal sphere) profile. Such 
conspicuous RDP irregularities suggest significant profile 
erosion or dispersion of stars from a primordial clust er (e.g. 
iBica. Bonatto fc Dutrall2008l ; IBonatto fc Bicall2009bl h 
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Figure 1. The optical environments are shown in 10' X 10' DSS B 
images of Pismis5 (top-left), vdB80 (right), NGC 1931 (bottom- 
left) and BDSB96 (right). Gas emission, dust reflection and/or 
absorption are present in all fields in varying proportions. Orien- 
tation: North to the top and East to the left. 
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Figure 2. 4' X 4' 2MASS K s images showing the embedded clus- 
ters Pismis5 (top-left), vdB 80 (right), NGC 1931 (bottom-left) 
and BDSB96 (right). Right ascension and declination are given 
in the abscissa and ordinate, respectively. Orientation: North to 
the top and East to the left. 



In some cases, early star cluster dissolution may lead to 
the formation of low-mass OB associations, the subsequent 
disp ersion of which may be an important source of field stars 
(e.g. iMassev, Johnson fc Gioia-Eastwoodlll995l ). Bochuml, 
with an irregular and clumpy RDP, can be an example of 
such an evolving structure. In this context, NGC 2244 ap- 
pears to be another example of a young OC in the process of 
dissolving in a few 10 7 yr. Perhaps the difference between ob- 
jects like Bochum 1 and NGC 2244 and normal young OCs 
(i.e. with the combined population of MS and PMS stars 
distributed according to a cluster RDP as in NGC 6611 and 
NGC 4755) is related to the interplay between environment 
conditions, star-formation efficiency and stellar mass. In- 
terestingly, the mass stored in the MS+PMS members of 
Bochum 1 and NGC 2244 is a factor 2-3 lower than those of 
NGC 6611 and NGC 4755. 

In this paper we investigate the nature of the 
poorly-studied, young and low-mass OCs Pismis 5, vdB 80, 
NGC 1931 and BDSB 96, by means of their photometric 
and structural properties. Their location in the Galaxy, 
in the 2 nd and 3 rd quad rants (Table [D, minimises field- 
star contamination (e .g. iBica, Bonatto fc Camargcl 120081: 
iBonatto fc Bicall2008bl ), which is essential when PMS stars 
are expected to dominate in Colour-Magnitude Diagrams 
(CMDs). Our main goal is to determine whether such young 
and low-mass systems can be characterised as typical OCs 
or if they are on their way to dissolution. In addition, we will 
derive their fundamental and structural parameters, most of 
these for the first time. 

This paper is organised as follows. In Sect. [2] we recall 
literature data on the target objects, discuss the 2MASS 
photometry and build the field-star decontaminated CMDs. 
In Sect. [3] we derive fundamental cluster parameters. In 
Sect. [4] we derive structural parameters. In Sect. [5] we esti- 



mate cluster mass. In Sect. [6] we compare structural param- 
eters and dynamical state with those of a sample of nearby 
OCs. Concluding remarks are given in Sect. [7] 



2 THE OCS AND THE 2MASS PHOTOMETRY 

The optical environments of the objects are shown in 10' x 
10' DSS3 B images (Fig.[T}, in which nebular gas and/or dust 
emission is present in varying proportions. The embedded 
clusters show up in the 4' x 4' 2MASS K s images (Fig. [2]). 

Pismis 5 (also known as ESO 313-SC7) is located in Vela 
and its field contains a few relatively bright stars mixed with 
gas emission and du st (Fig .[ T] top -left panel). It was discov- 



ered by P. Pismis (|Pismisl 19591 ) on Schmidt plates from 



the Tonantzintla Observatory (Mexico). According to SIM- 
BAE0, the 2 brightest stars in this field are CD - 39 4582 
(a(2000) = 08 h 37 m 38.0 s , 5(2000) = -39°35'02", V = 9.9, 
J = 9.4, K s = 9.1) and CD - 39 4599 (a(2000) = 
OS^SS™!^ 8 , 5(2000) = -39°34'6.1", V = 10.0, J = 10.0, 
K s — 9.9). Both appear to be members of Pismis 5 (Fig. [5}. 
As far as we are aware, the onl y work on Pism i s 5 is the 
XJBV+Hp photometric survey bv lVogt fc Moffatl l|l973l ). in 
which it was not considered as a star cluster. 

Cradled inside the reflection nebula vdB-RN80 
dvan den Berghl Il96rj ) in Monoceros, the field of the star 
cluster vdB 80 (Fig. [1] top-right panel) is similar to that 
of Pismis 5. The brightest star NSV 2998 (SIMBAD: 
a(2000) = 06 ,l 30 m 49.8 a , 5(2000) = -09°39'14.8", V = 8.8, 
J = 8.1, K s = 8.1) appears to be the only MS member of 



2 Extracted from the Canadian Astronomy Data Centre (CADC) 
- http://cadcwww.dao.nrc.ca/ 
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Table 1. Fundamental parameters of the target OCs 
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vdB 80 (Fig. P. In the only reference to vdB 80 as a star clus- 
ter, [Ahumada^elaL] l|200ll ) estimated the age 4.5 ± 1.5 Myr 
by means of integrated spectra. 

The star cluster NGC 1931 (Stock 9, Collinder68) lies 
within the bright nebula NGC 1931 (Sh 2-237) in Auriga. 
iMoffat. Fitzgerald fc Jackson! (Il979l ) determined the dis- 
tance from the Sun cLq = 1.8 kpc, while iPandev fc Mahral 
1 19861 ) found d Q = 2.2 kpc. iBhatt et all l| 19941 ) found dp, = 
2.2 kp c, E(B - V) = 0.55 and the age ~ lOMvr. lLata et all 
(|2002|) derived E(B - V) = 0.58 and ~ 13 Myr. Finally, 
IChen. Chen fc Shul (|2004l ) found d e = 3.1 kpc, ~ 10 Myr 
and the cluster radius R = 4'. 13. The latter work was based 
on 2MASS photometry, as in the present study with our 
tools, since the cluster is embedded in nebular emission. 

Discovered in a survey of infrared embedde d star clus- 
ters a nd stellar groups carried with 2MASS by IBica et all 
(2003), BDSB 96 (in the nebula Cederblad90, also known 
as Gum 3, Sh 2-297, RCWla, vdB-RN94 and Ber 134) in 
Canis Major, was classified as an infrared OC located at 
about 1.1 kpc from the Sun. Its field contains the single 
bright star HD 53623 (SIMBAD: a(2000) = 07 h 05 m 16.7 s , 
(5(2000) = -12°19'34.5", V = 8.0, J = 8.0, K s = 8.0), which 
is in the MS of BDSB 96 (Fig.0. Nebular gas and/or dust 
emission is more conspicuous in its field (Fig.[T] bottom-right 
panel) than in Pismis5, vdB 80 and NGC 1931. 

Table [1] provides parameters found in the literature and 
derived here. Since accurate spectral types of the bright stars 
listed above are not available, spectroscopic distances cannot 
be computed. The central coordinates were recomputed to 
match the absolute maximum present in the stellar surface 
densities (Sect. [2]3}- 



2.1 2MASS photometry 

The present OCs still retain part of the primordial gas and 
dust (Fig.[T]), which makes the near-IR the optimal spectral 
range to probe them. For instance, the number of detected 
stars at a given radius in NGC 2244 in the optical is signifi- 
cantly lower t han in the near-IR, dro pping to about 5% for 
the full field (jBonatto fc Bicall2009bh . For this purpose we 
work with the 2MASS3 near-IR J, H, and K 3 photometry, 
which allows the spatial and photometric uniformity useful 
for wide extractions that, in turn, provide high star-count 
statistics. Within this perspective, we have been developing 



quantitative tools to statistically disentangle cluster evolu- 
tionary sequences from field stars in CMDs. Decontaminated 
CMDs, in turn, have been used to investigate the nature of 
cluster c andidates and derive their astro physical parame- 
ters fe.g. IBica. Bonatto fc Camargoll2008l ). In short, we ap- 
ply (i) field-star decontamination to uncover the intrinsic 
CMD morphology, essential for a proper derivation of red- 
dening, age, and distance from the Sun, and (ii) colour- 
magnitude filters, which are essential for intrinsic stellar 
RDPs, as well as luminosity and mass functions (MFs). In 
particular, the use of field-star decontamination in the con- 
struction of CMDs has proved to constrain age and distance 
more than working with ra w (observed) photom etry, espe- 
cially for low-latitude OCs (jBonatto et al.|[2006al ). 

2MASS can reach adequate CMD depths for nearby 
OCs. For instance, our group has studied the young OCs 
NGC 6611, NGC 4755, NGC 2239 and NGC 2244. Abundant 
PMS stars were seen in the w 1 Myr old NGC 6611 (which 
is essentially embedded) and the 1 - 6 Myr old NGC 2244, 
and a few remaining ones in the « 14 Myr old NGC 4755. As 
nearby old er OCs we cite NGC 2 477 l|Bonatto fc Bicall2005l ) 
and M67 (|Bonatto fc Bicall2003l ). 

Photometry was extracted in a wide circular field of 
radius R cxt (Table [1} with VizieBQ. The extraction radii 
are large enough to allow determination of the background 
level (Sect. |3} and to statistically characterise the colour and 
magnitude distribution of the field stars (Sect. [!E2"jl . As pho- 
tometric quality constraint, only stars with J, H, and K s er- 
rors lower than 0.1 mag were considered- Reddening correc- 
tions are based on the relations Aj/Ay = 0.276, Ah/Av = 
0.176, AkJAv = 0.118 and A j = 2.76 x E(J - H) given by 
iDutra. Santiago fc Bical (12002)1. with Ry = 3.1, considerin g 
the extinction curve of ICardelli. Clayton fc Mathisl (|1989l ). 



2.2 Field decontamination 

Field-star decontamination is a very important step in the 
identification and characterisation of star clusters, espe- 
cially the clusters near th e Galactic equat o r. Diff erent ap- 
proaches are described in iBonatto fc Bical l|2009bl ). In this 
paper we employ the de c ontamination algorithm det a iled in 
IBonatto fc Bical J2007al). [Bica. Bonatto fc Camargol |2008) 



and lBica fc Bonattd 1 20081 ). and briefly described below. 



4 The Two Micron All Sky Su rvey, 
Sky data release (Skrutskic et al.1 Il997l) 

1:/ '/www. ipac.caltech. edu/2mass /releases/ allsky/ 



5 http://vizier.u-strasbg.fr/viz-bin/VizieR?-source= 11/ 24 6 
All 6 In all cases, the bright stars have photometric errors within the 

adopted range. Besides, any bias against low-mass stars (larger 
errors) introduced by this criterion (applied equally for cluster 
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2.3 Decontaminated surface density maps 





Figure 3. Stellar surface-density <r(stars arcmin -2 ) computed 
with field-decontaminated photometry to enhance the clus- 
ter/background contrast. A(o cos(<5)) and AS in arcmin. 



Figure [3] shows the spatial distribution of the stellar surface- 
density (a, in units of stars arcmin -2 ) around the 4 objects, 
measured with 2MASS photometry. We use field-star de- 
contaminated photometry (Sect. \Tl2^ to maximise the clus- 
ter/background contrast. The surface density is computed 
in a rectangular mesh with cells 2.5' x 2.5' wide, reaching 
total offsets of |A(a cos(5))| = |A<5| « 40' with respect to 
the cluster centre (TableQJ. Despite the gas and dust associ- 
ated with the present low-mass clusters (Fig.[l|, the central 
excesses are conspicuously detected in the decontaminated 
surface-density distributions. Besides, the residual surface- 
density around the centre has been reduced to a minimum 
level. 

By design, our decontamination depends essentially on 
the colour-magnitude distribution of stars located in dif- 
ferent spatial regions. The fact that the decontaminated 
surface-density presents a conspicuous excess only at the as- 
sumed cluster position implies significant differences among 
this region and the comparison field, both in terms of colour- 
magnitude and star counts within the corresponding colour- 
magnitude bins. This meets cluster expectations, which can 
be characterised by a single-stellar population, projected 
against a Galactic stellar field. 

The respective isopleth surface maps are shown in 
Fig- HI m which cluster size and geometry can be observed. In 
all cases, the central region (R < 5') appears essentially cir- 
cular, with elongated external regions, especially NGC 193f . 



Figure 4. Similar to Fig. [3] for the isopleth curves, for a better 
visualisation of cluster extent and geometry. 



The algorithm (i) divides the whole range of CMD mag- 
nitude and colours into a 3D grid of cells with axes along the 
J magnitude and the (J — H) and (J — K s ) colours, (ii) es- 
timates the number density of field stars in each cell based on 
the number of comparison field stars with similar magnitude 
and colours as those in the cell, and (Hi) subtracts the esti- 
mated number of field stars from each cell. Typical cell di- 
mensions are AJ = 1.0 and A( J - H) = A( J - K s ) = 0.2; 
the comparison fields are located within R = 30' and the ex- 
traction radius. With this setup, the subtraction efficiency, 
i.e. the difference between the background contamination 
(which may be fracti onal) and the number of subtracted 
stars in each cell (e.g. jBonatto fcF iica 2008b), summed over 
all cells, is higher than 90% in all cases. 



3 FUNDAMENTAL PARAMETERS 

J x (J — K s ) CMDs built with the raw photometry of the 
sample objects are shown in the top panels of Figs. [5] and 
[6] In all cases, the sampled region contains most of the 
cluster stars (Fig. 1 10 p . When qualitatively compared with 
the CMDs extracted from the equal-area comparison 
(middle panels), features typical of very young OCs emerge: 
a relatively vertical and poorly-populated MS together with 
a large number of faint and red PMS stars. 

The decontaminated CMDs are shown in the bottom 
panels of Figs. [S] and [5] As expected, essentially all contami- 
nation is removed, leaving stellar sequences typical of mildly 
reddened, young and low-mass OCs, with a developing MS 
and a significant population of PMS stars. 

Colour distributions wider than the spread predicted 
by PMS models occur in all cases (Fig. [jj, which reflects 
differential reddening. To examine this issue we show in 
Figs. HE reddening vectors computed with the 2MASS ra- 
tios (Sect. [2~T)l for visual absorptions Av = to 5. Taking 
into account the PMS isochrone fit (Fig. [7]), the differential 
reddening appears to be lower than Av = 5. 

We base the fundamental parameter derivation on the 
field-decontaminated CMD morphologies (Fig. [7J , using as 
constraint the combined MS and PMS star distribution. We 
adopt solar metallicity isochrones because the clusters are 
young and located not far from the Solar circle (see below), 



and field stars) is probably minimised by the field decontamina- 
tion (Sect. HOI. 



7 The equal-area field extractions serve only for qualitative com- 
parisons, since the decontamination uses a large surrounding area 



(Sect. [H). 
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Figure 5. 2MASS CMDs of Pismis5 (left) and vdB80 (right). 
Top panels: observed photometry. Middle: equal-area compari- 
son field CMDs. Bottom: decontaminated CMDs. Bright stars 
in common with SIMBAD are: CD - 39 4582 (large circle) and 
CD -39 4599 (diamond) in Pismis5, and NSV 2998 (large circle) 
in vdB 80. Reddening vectors for Ay = — 5 are shown in the 
bottom panels. 



a reg ion essentially occupied by [Fe/H] ~ 0.0 OCs l)Friell 
Il995l i. 

To represent the MS we use Padova isochrones 
|Girardi et alj|2002h computed with the 2MASS J, H, and 
The similar decontaminated CMD morphologies, 
typical of young ages (Fig. [7}, indicate similar age-solutions 
for the present objects. 

Sophisticated CMD fit approa ches are available, espe - 
cially for the MS (as summarised in lNavlor fc Jeff ries 2006). 
However, given the poorly-populated MSs, the 2MASS pho- 
tometric uncertainties for the lower sequences and the im- 
portant population of PMS stars, we decided for the direct 
comparison of isochrones with the decontaminated CMD 
morphology. Thus, fits are made by eye, taking the combined 
MS and PMS stellar distributions as constraint, allowing as 
well for differ ential reddening and photometric uncertainties. 
Isochrones of ISiess. Dufour fc Forestinil (2000) with ages in 
the range 0.2 — 20Myr are used to characterise the PMS se- 
quences. The results are shown in Fig.[7]and discussed below. 
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Figure 6. Similar to Fig. [5] for NGC 1931 (left) and BDSB 96 
(right). The bright star in BDSB 96 is HD 53623. 



Pismis5: Given the poorly-populated MS, acceptable 
fits to the decontaminated MS morphology are obtained 
with any isochrone with age in the range 1 — 10 Myr. The 
PMS stars are basically contained within the 0.2 Myr and 
10 Myr isochrones as well, which implies a similar age range 
as the MS. Thus, we take the 5 Myr isochrone as representa- 
tive solution, with a 4 Myr spread. In this case, CD — 39 4599 
is in the MS while CD — 39 4582 would be a massive 
(~ 5M Q ) PMS almost reaching the MS. 

With the adopted solution, the fundamental param- 
eters of Pismis5 are the near-IR reddening E(J — H) = 
0.13 ± 0.01 (E(B - V) = 0.42 ± 0.03 or A v = 1.3 ± 0.1), 
the observed and absolute distance moduli (m — M)j = 
10.4 ± 0.1 and (m - M)o = 10.04 ± 0.10, respectively, and 
the distance from t he Sun dp, = 1. ± 0.1 kpc. We adopt 
Rq = 7.2 ± 0.3 kpc (|Bica et al.ll2006t ) as the Sun's distance 
to the Galactic centre to compute Galactocentric distance^. 
For Rq = 7.2 kpc, the Galactocentric distance of Pismis 5 is 
Rgc = 7.5 ± 0.1 kpc, which puts it w 0.3 kpc outside the 
Solar circle. This solution is shown in Fig. [7] 

vdB 80: The bright star NSV 2998 (J=8.153) helps 
constrain the MS age t o 5 ± 2 Myr, which a grees with 
the 4.5 ± 1.5 Myr age of lAhumada et al l (|200ll ). Similarly 
to Pismis 5, the PMS stars are essentially found within 



8 http://stev.oapd.inaf.it/cgi-bin/cmd These isochrones 
arc very similar to the Johnson-Kron-Cousins ones (e.g. 
[Bessel_& i F5rctt 1988), with diff erences of at most 0.01 in colour 
llBonatto. Bica fc Girardill2004l) . 



9 Derived by mean s of the Glob u lar C luster (GC) spatial dis- 
tribution. Rec ently. iTrippe et all d2008l) found R GC = 8.07 ± 
0.32 kpc while iGhez et al.l fcoolj) found i? GC = 8.0 ± 0.6 kpc or 
Rgc = 8.4 ± 0.4 kpc, under different assumptions. 
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Figure 7. Adopted MS+PMS isochrone solutions to the decon- 
taminated CMDs. SIMBAD bright stars and reddening vectors 
as in Figs. [5] and [6] Shaded polygons show the MS (dark-gray) 
and PMS (light-gray) colour-magnitude niters (Sect.|?]l. 
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Figure 8. Decontaminated CMD of vdB80, used to illustrate 
the PMS mass estimate. Light lines: MS and PMS isochrones set 
as in Fig. [7] Heavy-dotted: PMS evolutionary tracks for different 



the 0.2 — lOMyr isochrones. With this solution, the funda- 
mental parameters of vdB80 are E(J - H) = 0.19 ± 0.03 
(E(B - V) = 0.61 ± 0.10 or A v = 2.0 ± 0.3), (to - M)j = 
12.1 ± 0.3, (m - M)o = 11.58 ± 0.31, d Q = 2.1 ± 0.3 kpc, 
and Rgc = 8.9 ± 0.3 kpc, thus m 1.7kpc outside the Solar 
circle. 

NGC 1931 : Besides the most populated MS of the present 
sample, NGC 1931 presents the smallest gap between the 
faintest MS star and the PMS distribution, which suggest 
a more evolved phase. Indeed, the PMS stars distribute be- 
tween the 0.2 — 20 Myr isochrones. Based on this, the MS age 
ca n be set at 10 ± 3 Myr, which agrees with t h e estimates 
oflPandev fc Mahral I dl986l). iBhatt et alj (| 19941 ). [Lata et all 
|2002l ) and lChen. Chen fc Shul (|2004l ) . This solution implies 
E(J -H) = 0.19 ± 0.03 (E(B - V) = 0.61 ± 0.10 or A v = 
2.0±0.3), (m - M)j = 12.4±0.3, (m - M)o = 11.86±0.22, 
d© = 2.4 ± 0.2 kpc, and R GC = 9.6 ± 0.2 kpc, thus « 2.4 kpc 
beyond the Solar circ le. Our value for dp, is intermediate 



between the 1. 8 kpc oflMoffat. Fitzgerald fc Jackson! ijl979h 
and 3.1 kpc of lChen. Chen fc S hu (2004J), and agrees with 
the 2 .2 kpc of iPandev fc Mahral (|l986l ) and IBhatt et all 
|l994h . 

BDSB 96 : The CMD of this object is similar to that of 
vdB 80 - including the presence of a bright star (HD 53623) 
in the MS, which suggests a similar solution, with the MS 
age 5±3Myr. This solution implies E{J - H) = 0.12 ±0.03 
(E(B - V) = 0.37 ± 0.10 or A v = 1.1 ± 0.3), (to - M)j = 



11.0 ± 0.2, (m - M) = 10.68 ± 0.22, d Q = 1.4 ± 0.2 kpc, 
and Rgc = 8.2 ± 0.2 kpc, thus ~ 1.0 kpc outsid e the Solar 
circle. Our value of ci© agrees with the 1.1 kpc of iBica et al.l 
(|2003t ). 

Besides the MS age, the present objects have in common 
a significant age spread implied by the PMS stars. This indi- 
cates a non-instantaneous star formation process, similar to 
what was found in our previous studies of, e.g. NGC 4755, 
NGC 6611 and NGC 2244. 



3.1 Colour-colour diagrams 

When transposed to the near-IR colour-colour diagram 
(J — K s ) x (H — K s ) (Fig. |9}, the age and reddening so- 
lutions derived for the present objects consistently match 
their field-star decontaminated p hotometry. Since they in- 
clude PMS stars, we use tracks of[g iess. Dufour fc Fore stini 
(2000) to characterise the age. MS stars lie on the blue 
side of the diagram. As expected from the CMDs, a sig- 
nificant fraction of the stars appears to be very reddened. 
Besides, since most stars have (H — K s ) colours close to the 
isochrone, within the uncertainties, the fraction of stars still 
bear ing circunstellar discs (with excess in H) appears to be 
low l|Bonatto et al.|[200 6b and references therein). 
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Figure 9. Decontaminated colour-colour diagrams. PMS 
isochrones set with the derived reddening values (Sect. SIM- 
BAD bright stars and reddening vectors as in Figs. [5] and [6] 



4 CLUSTER STRUCTURE 

We use the RDPs, defined as the projected stellar number 
density around the cluster centre (i.e. the maximum of the 
surface density maps - Sect. l2.3] l. to derive structural param- 
eters. To minimise noise, we work with colour-magnitude 
filters (shown in Fig. [7} to exclude stars with colours un- 
like those of the cluster CMD morpholog-J^l. This enhances 
the RDP contrast relative to the b ackground, especially in 
crowded fields (e.g. lBonatto fc Bicall2007al ) . Examples of the 
advant age of using colour-ma g nitude filters can be found 
in, e.g. lBonatto fc Bica (l2007bT). Bonatto. Bica fc Santos Jr.l 



(|2008h and IBonatto fc Bical (|2008bT r 

Rings of increasing width with distance from the cluster 
centre are used to preserve spatial resolution near the cen- 
tre and minimise noise at large radii. The set of ring widths 
used is AR = 0.25, 0.5, 1.0, 2.0, and 5', respectively for 
0' «: R < 0.5', 0.5' ^ R < 2', 2' R < 5', 5' < R < 20', and 
R > 20'. Because of the low number of stars in the central 
parts of vdB 80, we used AR = 0.5 for R < 1'. The residual 
background level of each RDP corresponds to the average 
number-density of field stars. The R coordinate (and uncer- 
tainty) of each ring corresponds to the average position and 
standard deviation of the stars inside the ring. As a caveat 
we note that the present OCs are not spherical, especially 



at the outskirts (Fig. [4} , which might affect the RDPs for 
large radii. Because of this, deviations in the central parts 
of the RDP are not expected to be significant. 

Minimisation of the number of non-cluster stars by the 
colour-magnitude filters yielded RDPs (Fig. I10|l highly con- 
trasted relative to the background. For simplicity we fit the 
RDPs with a(R) = a bg + cr /(l + {R/R c ) 2 ), where a bg is 
the residual background density, cro is the central density of 
stars, and R c is the core radius. This func tion, applied to 
star counts, is similar to that introduced by iKind (|l962l ) to 
describe the surface-brightness profiles in the central parts 
of gc£B To minimise degrees of freedom, cro and R c fol- 
low from the fit, while Ob g is measured in the field and kept 
fixed. The best-fit solutions and uncertainties are shown in 
Fig. 1101 and the parameters in Table [2] We also estimate 
the cluster radius (-Rrdp), i-e. the distance from the centre 
where the cluster RDP an d field fluctuations ar e statisti- 
cally indistinguishable (e.g. iBonatto fc Bical l2005T l. and the 
density contrast parameter 8 C — 1 + ao/ubg (Table [2}. Inter- 
estingly, the density-contrast parameter reaches high values, 
7 ^ 23, as expected from compact star clusters. 

Within uncertainties, the adopted King-like function 
describes the RDPs of vdB 80 and BDSB 96 along the whole 
radius range. To a lesser degree, the same applies for Pismis 5 
and NGC 1931, in which the innermost bin presents an ex- 
cess over the fit. In old star clusters, such a cusp has been at- 
tributed to a post-c ore collapse structure, like those d etected 
m some GCs (e.g. iTrager. King fc Diorgpvskilfl995n . Gyr- 
old OCs, e.g. NGC 3960 l|Bonatto fc Bicall2006n and LK 10 
l|Bonatto fc Bicall2009al ). also display such feature, which is 
related to dynamical evolution. Wi th respect to very young 
clusters, the R DPs of NGC 2244 feonatto fc Bical l2009bf ) 
and NGC 6823 l|Bica. Bonatto fc DutralbOOct ) alsopresent a 
central cusp, similarly to those in Pismis 5 and NGC 1931. 
Clusters are not expected to dynamically evolve into a post- 
core collapse on such short time-scales. Consequently, the 
cusp in young clusters must be related to molecular cloud 
fragmentation and/or star- forming effects, and may suggest 
early deviation from dynamical equilibrium (Sect.|SJ). 

Taken at face value, the core radii of the present ob- 
jects (0.2 < i?c(pc) < 0.5) fall on the low-7? c tail of the 
dis tribution derived for a sample of relatively nearby OCs 
bv lPiskunov et al l (|2007l ). 



5 MASS ESTIMATE 

The most conspicuous CMD feature of the present OCs is 
the dominant presence of PMS stars, followed by the poorly- 
populated MS (Fig. 0. Obviously, in each cluster, most of 
the mass is stored in the PMS stars. Because the decontam- 
ination algorithm actually excludes stars (in integer num- 
bers), the decontaminated photometry should be used es- 
sentially to determine the CMD morphology or the colour- 
colour distribution (Sect. [3}. However, when magnitude (or 
alternatively mass) bins are used to build the MFs, the bin- 
to-bin subtraction of the comparison field (normalised to 



10 They are wide enough to include cluster MS and PMS stars, 
together with the photometric uncertainties and binaries (and 
other multiple systems). 



11 Besides, because of the relatively small number of stars in 
the present OCs, fluctuations in surface-brightness profiles are 
expected to be larger tha n those in RDPs . Alter native RDP fit 
functions are discussed in IBonatto fc Bical ll2008al) . 
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Table 2. Derived structural parameters 



Cluster 




TO 


Rc 


Rrdp 


Sc 


1' 




CO 




-Rrdp 




(*'~ 2 ) 


(*'- 2 ) 


(') 


(') 




(pc) 


(*pc~ 2 ) 


(*pc~ 2 ) 


(PC) 


(PC) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Pismis 5 


1.81 ±0.02 


10.2 ±4.3 


0.69 ±0.22 


6.0 ±0.3 


6.6 ± 2.3 


0.296 


20.7 ±0.3 


116.7 ±49.2 


0.20 ±0.06 


1.8 ± 0.1 


vdB 80 


1.25 ±0.02 


20.3 ± 5.0 


0.46 ± 0.08 


5.8 ± 0.3 


17.2 ±4.0 


0.599 


3.5 ±0.1 


56.4 ± 14.0 


0.28 ±0.05 


3.5 ± 0.2 


NGC 1931 


2.08 ± 0.02 


19.6 ± 6.1 


0.70 ±0.15 


8.5 ± 0.5 


10.4 ± 2.9 


0.684 


4.4 ±0.1 


41.9 ± 13.0 


0.48 ±0.10 


5.8 ± 0.3 


BDSB 96 


1.29 ±0.01 


28.2 ± 11.8 


0.88 ± 0.30 


11.0 ± 1.0 


22.9 ±9.2 


0.397 


8.1 ±0.1 


178.8 ± 74.8 


0.35 ±0.12 


4.4 ±0.4 



Table Notes. Core (R c ) and cluster (Rrtjp) radii are given in angular and absolute units. Col. 6: cluster/background density contrast 
parameter (<5 C = 1 + ao/a^g). Col. 7: arcmin to parsec scale. 



10 



10 



Pismis5 

T 

+ 


vdB80 




R »" T\ 




\ R RDP 






: — — — 


: — — — 


I NGC1931 


j BDSB96 


T 






\ R RDP 















0.1 



1 10 0.1 

R (arcmin) 



1 10 

R (arcmin) 



Figure 10. Stellar RDPs built with colour-magnitude filtered 
photometry together with the best-fit King-like profile (solid line), 
the lcr uncertainty (light-shaded region) and the background level 
(shaded polygon). Note the central cusp in the RDPs of Pismis 5 
and NGC 1931. 



the same projected area as the cluster) contribution is ex- 
pected to produce fractional n umbers, which should be in- 
corporated into the MFs (e.g. iBonatto fc Bical I2005T ). For 
this reason we use the colour-magnitude filtered photome- 
try (Fig. [7} to estimate the cluster mass. 

The number of MS stars is derived by counting the stars 
within R ^ -Rrdp (in bins of A J = lmag), and subtracting 
those in the field (normalised to the same area). The corre- 
sponding stellar mass in each magnitude bin is taken from 
the mass-luminosity relation derived from the isochrone fit 
(Sect. [3]). Summing the values in each magnitude bin pro- 
duces the total number (ums) and mass (mMs) of MS stars. 



A similar strategy is applied to the PMS stars. We consider 
the evolutionary tracks for the PMS masses 0.5, 1, 2, 3, 5 
and 7 M Q (Fig. [£} . The number of PMS stars between any 
2 tracks is counted for R ^ Rrdp and the field, taking into 
account the reddening vectors. The average mass between 
two evolutionary tracks is taken as the mass of a respective 
PMS star. All tracks summed result in the number (upms) 
and mass (rnpMs) of PMS stars (Table [3}. 

As anticipated by the CMDs (Fig. 0, the MSs as a 
rule, are poorly-populated. Pismis 5, vdB80 and BDSB 96 
contain about 4 MS stars, while NGC 1931 hosts about 
15 MS stars. Such low numbers reflect the young age and 
low-mass nature of these objects. Indeed, the mass of the 
present objects, as derived from their MS±PMS mem- 
bers, are within 60- 180MeI3. For comparison, Bochuml, 
NGC 6823, NGC 2244 and NGC 2239 contain 128, 65, 26 and 
26 MS stars, respectively, while their (MS±PMS) masses are 
720, 1150, 625 and 301 M , respectively. 

An estimate of the MS and PMS MFs can be made 
with the approach described above. In Fig. [11] we show the 
resulting MFs for the PMS and MS stars combined. Since 
the MSs are very-poorly populated (Table [3J , it should be 
noted that these MFs correspond essentially to the PMS 
stars; the MSs contribute somewhat to the range m > 2 Mq. 
Besides, given the simplifying assumptions, the error bars 
in the MFs are only formal, and should be taken as a lower 
limit. Bearing this in mind, the MFs can be represented by 
the function 4>{m) = g oc m" (1+x) , wit h slopes ^ (Tab led 



significantly flatter than the % ~ 1-35 of ISalpeterl l l 19551 ) ini- 
tial mass function (IMF). Interestingly, these MFs do not 
appear to present a turnover at the sub-Sol ar m ass range, 
as sug gested by the IMFs of iKroupal l|200ll ) and IChabrierl 
(2003). However, given the low number of stars and the 
simplifying assumptions, we have omitted several possible 
systematic biases associated with the IMF determination: 

(i) uncertainties in the low- mass IMF due to differential 
reddening and possibly crowding (the latter should be min- 
imal, since the objects are intri nsically poorly popul ated), 

(ii) missing c ompanion stars (e.g. Maiz Apella niz 2008]), (in) 



s ( e.g. | Ma 

binning (e.g. iMaiz Apellanidl2005l ). and (iv) errors on the 
mass/luminosity relation of individual stars. 



12 As a caveat we note that, because of the differential reddening 
and the 2MASS photometric limit, we may not detect the PMS 
stars of very-low mass; thus the mass values may be somewhat 
higher. 
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Table 3. MS and PMS stellar content 







MS 




PMS 




MS+PMS 




Cluster 


A m 


nMS 


m M s 




mpAfs 


riMS+PMS 


mus+PMS 


X 




(M e ) 


(stars) 


(Mq) 


(stars) 


(M Q ) 


(stars) 


(Mo) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Pismis 5 


1.9-13.4 


3±1 


12 ±3 


100 ± 20 


46 ±7 


103 ± 21 


58 ± 8 


-0.03 ±0.19 


vdB 80 


2.3-21.6 


4± 1 


40 ±11 


108 ± 12 


55 ± 13 


112 ±13 


95 ± 17 


-0.24 ±0.09 


NGC 1931 


2.0-17.8 


15 ±3 


62 ± 13 


183 ± 33 


115 ±37 


198 ± 34 


177 ± 39 


-0.02 ±0.08 


BDSB 96 


1.8-23.8 


5±1 


23 ±6 


171 ± 23 


131 ±21 


176 ± 24 


154 ± 23 


+0.80 ±0.08 



Table Notes. Col. 2: Effective MS mass range. Col. 9: MS+PMS mass function slope \, derived from the fit of 4>(m) tx m 
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Figure 11. MS+PMS mass functions fitted by the function 
<f>(m) oc m^( 1 + x ' (dashed line). 



6 DISCUSSION 

To put Pismis 5, vdB 80, NGC 1931 and BDSB 96 into per- 
spective, their astrophysical parameters (Sects. |3j [4] [5} 
are compared to those of a sample of nearby OCs in dif- 
ferent environments. The reference sample contains some 
bright nearby OCs (|Bonatto fc Bical 120051 ; iBonatto et all 
l2006bl ) together with a group of OCs projected towar ds the 
central parts of the Galaxy (Bonatto & Bica 2007a). The 
young OCs NGC 6611 llBonatto. Santos Jr~fc Bical l200r3 ). 
NGC 6823 jBica Bonatto fe Dutral 120081 ) and NGC 2239 
ijBonatto fc Bical l2009bi r are included for comparison 
with gravit ationally bound objects of si milar ages, while 
Bochum 1 jBica. Bonatto fe Dutral l2008h and NGC 2244 




10 10 10 10 10 
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Figure 12. Diagrams dealing with astrophysical parameters of 
OCs. Gray-shaded circles: reference OCs. The young star clusters 
NGC 6611, NGC 6823, Bochum 1, NGC2244 and NGC 2239 are 
indicated for comparison purposes. Analytical relations in panels 
(d) and (e) are discussed in the text. Dotted line in panel (c): 
Salpeter's IMF slope \ = 1-35. 



i Bonatto & Bica 2009b) might be dynamically evolving into 
OB associations. All objects have been similarly analysed. 

We work with the diagnostic-diagrams shown in Fig. 1121 
Panels (a) and (b) examine the dependence of cluster 
(Rrdp) and core (R c ) radii on age, respectively. Pismis 5, 
vdB 80 and BDSB 96 (and to a lesser degree NGC 1931) 
present small R c and .Rrdp, similar to NGC 6823 and typ- 
ical of Gyr -old OCs undergoing disruption near the Solar 
circle (e.g. iBergond. Leon fc Guilbertl l200ll ; lLamers et al.l 
2005). Within uncertainties, the 4 objects follow the rela- 
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tion R RD p = (8.9 ±0.3) x R^- 0±0 A) (panel d), derived with 
the reference sample. They also appear to follow the relation 
of increasing OC size with Galactocentric distancj^l (panel 
f)- 

When the mass-density r adial distributi o n fol- 
lows a King-li ke profile (e.g. Bonatto fe Bical 2007bl; 

Bonatto fe Bical l2008al ; iBonatto. Bica fe Santos Jr.l 



2008), the mass inside Rrdp is a function of 
R c and the central mass-surface density {(J mo), 
M c i us = TrRco-Moln [l + (-Rrdp/-Rc) 2 ]. With the relation 
between R c and -Rrdp (panel d), this equation becomes 
Mdus w 13.8(Jmo Rc- The reference OCs, together with the 
present sample, are contained within King-like distributions 
with c mo constrained within 30 < ctmo (M©pc~ 2 ) < 600 
(panel e) . The exception is NGC 2244, which appears to 
present too big a core for the total mass. 

Finally, when the total (MS+PMS) MF slope is con- 
sidered (panel c), Pismis5, vdB80 and NGC 1931 present 
MFs significantly flatter than those of similarly young OCs. 
To a lesser degree, the same applies to the newly formed 
BDSB 96. Such fiat slopes are equivalent to those derived for 
some very old OCs in the reference sample that, in general, 
unde rgo advanced dynamical evolution (e.g- lBonatto fe Bical 
120051 ). Given the young age, low mass, irregular RDP and 
the flat mass function, the OCs NGC 1931 and Pismis 5 may 
be evolving into OB associations or remnants in a few 10 7 yr. 
Since both OCs are young, such effects may be related to 
another mechanism than age-related dynamical evolution. 
It is probably associated with star formation. vdB 80 and 
BDSB 96, on the other hand, appear to be typical young 
and low-mass OCs. 



7 SUMMARY AND CONCLUSIONS 

Probably because of the interplay between environmental 
conditions, star-formation efficiency and stellar mass, only 
a few percent of the embedded clusters survive the first few 
tens of Myr. Thus, the derivation of astrophysical parame- 
ters of OCs in this phase may shed some light on the roles 
played by the above process in the dissolution/survival is- 
sue. In this context, poorly-populated and massive OCs in 
different environments are yet to be investigated. 

In the present paper we derive astrophysical parameters 
and investigate the nature of the young and low-mass OCs 
Pismis 5, vdB 80, NGC 1931 and BDSB 96. The wide-field 
and near-IR depth provided by 2MASS (with errors lower 
than 0.1 mag) are employed coupled to field-star decontami- 
nated photometry. This enhances cluster CMD evolution- 
ary sequences and stellar radial density profiles, yielding 
more constrained fundamental and structural parameters. 
The OCs are located within 174° < I < 260° and |6| < 9°, 
thus, the errors potentially induced by the background star 
subtraction are not critical. 

The decontaminated CMDs exhibit similar properties, 
basically a poorly-populated MS, a dominant fraction of 
PMS stars together with some differential reddening. MS 
ages are constrained within 5 ± 4 Myr (Pismis 5, vdB 80, 

13 It may be partly primordial, in the sense that the high molec- 
ular gas density in central Galactic regions may have pro duced 
small clusters (e.g. Ivan den Bergh. Morbev fe Paz dcr 1991). 



BDSB 96) and 10±3Myr (NGC 1931). However, the PMS 
stars suggest a wider age spread (~ 20 Myr), consistent 
with a non-instantaneous star formation process. The total 
(MS+PMS) stellar masses are low, within ~ 60 - 180 M©, 
with mass functions significantly flatter than Salpeter's 
IMF. 

The present OCs are rather small (R c < 0.48 pc, 
Rrdp J$ 5.8 pc), particularly Pismis 5 with R c ~ 0.2 pc and 
Rrdp ~ 1.8 pc. Structurally, the (MS+PMS) stellar RDPs 
follow a cluster-like profile for most of the radius range. Ex- 
ceptions are the inner regions of NGC 1931 and especially 
Pismis 5, which present a marked stellar-density excess. At 
~ 10 Myr, such a central cusp cannot result from large- 
scale cluster dynamical evolution. Instead, it probably is 
associated with molecular cloud fragmentation and/or star- 
formation effects. 

BDSB 96 and vdB 80 present structural properties of 
typical young, low-mass OCs, although with flat mass func- 
tions. On the other hand, the irregular RDPs - together with 
the low cluster mass and flat mass functions - of NGC 1931 
and especially Pismis 5, suggest that both OCs deviate from 
dynamical equilibrium. They are possibly evolving to be- 
come OB associations or remnants in a few 10 7 yr. 
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